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ABSTRACT 

We present high-angular resolution images of the post-AGB nebula 
IRAS 18276-1431 (also known as OH 17.7-2.0) obtained with the Keck II Adap- 
tive Optics (AO) system in its Natural Guide Star (NGS) mode in the Kp, Lp, 
and Ms near-infrared bands. We also present supporting optical F606W and 
F814W HST images as well as interferometric observations of the ^^CO (J=l- 
0), ^^CO (J=l-0), and 2.6mm continuum emission with OVRO. The envelope 
of IRAS 18276—1431 displays a clear bipolar morphology in our optical and 
NIR images with two lobes separated by a dark waist and surrounded by a 
faint 4'.'5x3'.'4 halo. Our i^'p-band image reveals two pairs of radial "search- 
light beams" emerging from the nebula center and several intersecting, arc-like 
features. From our CO data we derive a mass of M >0.38 [D/3 kpc]^ Mq and 
an expansion velocity \4xp=17km s~^ for the molecular envelope. The density 
in the halo follows a radial power-law (xr~^, which is consistent with a mass- 
loss rate increasing with time. Analysis of the NIR colors indicates the pres- 
ence of a compact central source of ~300-500 K dust illuminating the nebula 
in addition to the central star. Modeling of the thermal IR suggests a two- 
shell structure in the dust envelope: 1) an outer shell with inner and outer ra- 
dius _Rin~1.6xl0^^cm and -Rout ^1 -25 x 10^ ''cm, dust temperature Td~105-50K, 
and a mean mass-loss rate of M^10~^Mq yr~^; and 2) an inner shell with 
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i?in~6.3xl0^^cm, rdust~500-105K, and M^SxlO-^M© yr-^ An additional 
population of big dust grains (radius a>0.4mm) with Tdust =1 50-20 K and mass 
Mdust=(0. 16-1. 6)xlO~^[D/3kpc]^M0 can account for the observed sub-mm and 
mm flux excess. The mass of the envelope enclosed within i?out=l-25xl0^^cm 
derived from SED modeling is ~l[D/3kpc]^M0. 

Subject headings: instrumentation: adaptive optics, stars: AGB and post-AGB, 
stars: imaging, stars: mass loss, circumstellar matter, ISM: jets and outflows, 
planetary nebulae: general, infrared: stars 

1. Introduction 

In recent years, ~0'/l-resolution imaging with HST has provided many new insights into 
the evolution from the Asymptotic Giant Branch (AGB) to the phase of Planetary Nebula 
(PN). This transition, which takes place through the intermediate short-lived (~10^yr) stage 
of pre-planetary nebula (PPN), is characterized by dramatic changes of the nebular morphol- 
ogy and kinematics: the spherical, slowly expanding (V^xp'^^lSkm s~^) AGB circumstellar 
envelope (CSE) develops clear departures from sphericity and fast (> 100 km s~^) outflows 
directed in opposing directions along one or more axes. Although there is no consensus yet 
for what causes this spectacular metamorphosis, fast jet-like winds have been hypothesized 
to play an important role (see e.g. a recent review by Balick & Frank, 2002). In light of the 
many structural details revealed by HST data, e.g., the multiple lobes and the high degree 
of point symmetry of outflows observed in most young PNs, and PPNs, Sahai & Trauger 
(1998) hypothesized that fast coUimated outflows are the primary agent for setting the stage 
in the shaping of PNs, rather than a pre-existing equatorial density enhancement in the AGB 
CSEs, as assumed in the popular Generalized Interacting Stellar Winds model (Balick 1987). 
These outflows carve out a complex imprint within the AGB CSE producing and shaping 
the fast, bipolar lobes observed in most PPNs and PNs. 

In order to study the poorly known mass-loss processes mediating the AGB-to-PN 
transition, we have obtained 1.6-4.7;Lim images of PPNs with the Keck Adaptive Optics 
(AO) system (Sanchez Contrcras ct al., 2006). This forms part of our extensive, multi- 
wavelength studies of PPNs, which also include optical/NIR HST imaging and optical 
spectroscopy (Sahai & Sanchez Contreras 2004) and millimeter wavelength interferometry 
(Sanchez Contreras & Sahai 2004). Progress in our understanding of post-AGB evolution 
requires high spatial resolution because most PPNs are very small in angular extent {^5"), 
exhibit rich structure at a scale <0'/l, and because the shaping mechanism operates very 
close to the central star, <10^^cm (Balick & Prank 2002, Sanchez Contreras & Sahai 2001). 
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Also, observing at long wavelengths is desirable, since the dusty CSEs of most post-AGB 
objects are very optically thick. 

In this paper we present Kp, Lp, and AO images of the PPN IRAS 18276-1431, 
together with optical MVFPC2/HST and interferometric 2.6 mm (CO hne and continuum 
emission) data, which have helped our understanding of the complex structure of this object. 

Observational details and results are described in Sections 2 and 3. The analysis of the NIR 
and optical colors to derive the dust extinction towards the central star and the nature of the 
illuminating source are presented in Sections 4 and 5. Structural components and physical 
parameters (mass, mass-loss rate, and dust temperatiue) of IRAS 18276—1431 are discussed 
in Section 6 and 7. The formation and evolution of IRAS 18276—1431 is described in Section 
8. Our main conclusions are summarized in Section 9. 



1.1. The pre-planetary nebula IRAS 18276-1431 

IRAS 18276-1431 (also known as OH 17.7-2.0) is an OH/IR star that has reached the 
end of its evolution on the AGB and is rapidly evolving to the PN phase. IRAS 18276—1431 
displays the classical bimodal spectral energy distribution (SED) observed in PPNs, which 
is indicative of a detached CSE with relatively cool (~150K) dust grains (Le Bertre 1987). 
The progressive disappearance of H2O maser emission in IRAS 18276—1431 is consistent 
with a recent drop of the mass-loss rate from M>10~^ to well below M—10~'^Mq yr~^ 
(Engels 2002). The cessation of the large scale mass-loss, indicated by both the disappear- 
ance of the H2O masers and the presence of a detached envelope, is believed to designate the 
beginning of the post-AGB evolution. The non-detection of SiO masers in IRAS 18276—1431 
(Nyman et al. 1998) is also consistent with the post-AGB nature of this object. The central 
star of IRAS 18276—1431 has a spectral type earlier than ~K5, with an effective tempera- 
ture in the range reff=4000-10,000 K, and a total luminosity of L ~1.5xl0^ [D/A.6kpc]^ (Le 
Bertre et al. 1989). This is a confirmation that the object has left the AGB. 

Ground-based imaging polarimetry in the J and K near-infrared bands shows an ex- 
tended, non-spherical scattering envelope around IRAS 18276— 1431 (Gledhill 2005). The 
angular size of the envelope and position angle of its major axis in the J {K) band are 
~3"x2" and PA=24±2° (3'.'9x3'/4 and PA=22±2°), respectively. The polarized intensity 
images reveal a clear bipolar morphology and polarization vectors that are mostly aligned 
perpendicular to the bipolar axis, at PA^23°. This is typical of scattering in an optically 
thick envelope, where light escapes preferentially along bipolar cavities (Bains et al. 2003; 
Gledhill 2005). The axisymmetric scattering models by Gledhill (2005) suggest that the lobes 
of IRAS 18276—1431 lie very close to the plane of the sky, with an approximate inclination 



-4- 



of i~3°. 

Strong OH maser emission is observed in the equatorial regions IRAS 18276—1431 
(Bains al. 2003). This dense, equatorial component is elongated in the direction of PA~110°, 
i.e., is approximately perpendicular to the polar cavities, and has inner and outer angular 
radii of 0'.'25 and 0'.'85, respectively. As discussed by Bains et al. (2003), horn-hke fea- 
tures observed in the OH maser position-velocity diagram (see their Fig. 10) probably trace 
the base of a bipolar outflow very close to the central source. The relative absence of OH 
masers in the polar regions (compared with the much more intense OH emission observed 
in the equatorial regions) is consistent with the interpretation of these polar regions as cav- 
ities within which the lower gas density fails to generate strong maser emission. These OH 
maser observations also suggest a large-scale, regular magnetic field of strength B=-|-4.6mG 
(1612MHz) and B=+2.5mG (1667MHz) in IRAS 18276-1431 that could have played a role 
in causing the equatorially enhanced mass loss in this object. A burst of strongly polarized 
OH maser emission that could be linked to recent evolutionary changes in the shell has been 
reported by Szymczak & Gerard (2005). 

The distance to IRAS 18276—1431 is uncertain. OH phase-lag measurements indicate 
a distance in the range D=2-5.4kpc (Bowers et al. 1983; Herman & Habing 1985). The 
near kinematic distance is Z)=4.6kpc (Le Bertre et al. 1989). In this paper, we adopt an 
intermediate value of D—3 kpc. 

2. Observations 

2.1. Necir-Infrared Natural Guide Stcir AO imaging 

We have imaged the dust envelope around the OH/IR star IRAS 18276—1431 through 
the Kp(2.12/im), Lp(3.8;um) and M^(4.7/im) fihers with the W. M. Keck II 10-m telescope 
using the NIRC2 Camera (Matthews et al. 2005, in prep) with Adaptive Optics (AO) 
(Wizinowich et al. 2000). IRAS 18276-1431 was observed with AO in its Natural Guide 
Star (NGS) mode on UT 2005 August 12 at airmass 1.22-1.26. The AO guide star, USN- 
0755-0509985 (my=11.0) is located 23'/2 S and 0'/13 E from IRAS 18276-1431. The AO 
system ran at the speed of 315 Hz with an average intensity per sub-aperture of 120 counts. 
The field of view (FoV) and pixel scale of the NIRC2 narrow camera is 10" x 10" and 9.94 
mas pixel" ^. The seeing was 0'.'4 in the visible as measured on the acquisition camera. 

The images in the different bands were acquired using a generic 3-position dithering 
script avoiding the noisiest quadrant of NIRC2 where the number of frames for each position 
was 3 in Kp, 6 in Lp and 6 in Mg. The images have been reduced in a standard way. First, 
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we derived a sky frame by computing the filtered median of all frames for each band. Then, 
each frame was sky-subtracted, fiat-fielded, corrected for bad pixels and shift-and-added in 
a final image for each observing setup. We also observed the photometric standard #148 
S810-A from Lcggctt ct al. (2003) at UT 08:30 in the same region of the sky that we used to 
derive the photometry zero points for Kp, Lp, and Mg. We have used the magnitude-to-mJy 
flux density conversion factors given by Tokunaga and Vacca (2005). 

The spatial resolution in the images is nearly diffraction limited for the Keck 11 telescope 
in Kp (60-65 versus 48 mas) and diffraction-limited in Lp and Mg (Tabic 1). The PSF of our 
images has the typical elongation resulting from the AO anisoplanatism effect, as expected 
when the AO guide star is offset from the science target. This effect is seen in the Kp band 
images where the PSF varies from 55 mas near the AO guide star to 75 mas in the most 
distant imaged stars. This elongation is not seen in Lp and Mg as the isoplanatic angle 
becomes much larger (> 50"). We have used the nearest stars to IRAS 18276—1431 to check 
the image quality and characterize the PSF in our target (Table 1). The Strehl ratio (SR) 
was measured to be 0.16 (0.75) in the Kp (Lp) image, giving a scaled SR of 0.80 in Mg. 

2.2. Optical F606W and F814W HST imaging 

IRAS 18276-1431 was imaged on UT date 2002 Feb 16 (GO/SNAP program 9101) with 
the 800 X 800-pixel planetary camera of the Wide-Field & Planetary Camera 2 (WFPC2), 
which has a plate scale of O'.'0456/pixel, using the F606W [A = 0.60/im, AA = 0.123//m; 
exposures were 2x5 s and 2x100 s] & F814W [A = 0.80/i.m, AA = 0.149//m; exposures were 
2x3 & 2x60s] filters. For both filters, each of the two long-exposure images were dithered 
relative to each other by about 7.7 pixels to improve the image quality. The standard 
pipeline processed images were retrieved from the HST archive. The short-exposure images 
were used to check for saturated pixels in the long-exposure ones, and none were found. 
The long-exposure images were subsampled by a factor 2, and then registered to a common 
reference frame with high precision (i.e., to within <0.1 pixel) using field stars. Cosmic-ray 
removal and averaging was then carried out using the registered images. 

2.3. Millimeter-wavelength mapping with OVRO 

Interferometric mapping of the ^^CO J=l-0 and ^^CO J=l-0 emission hnes at 115 and 
110 GHz of IRAS 18276— 1431 was carried out using the millimeter interferometer of the 
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Owens Valley Radio Observatory (OVRO^) on Sept 29, 2002, with 5 antennas and baselines 
ranging from 18 m to 55 m. The units of the digital spectral line correlator were arranged 
to provide a total bandwidth of 90 MHz (~ 234 km s~^) with a channel spacing of 1 MHz 
(corresponding to ~2.6km s^^). The 3mm continuum emission from IRAS 18276—1431 was 
observed simultaneously using the dual-channel analog continuum correlator, which provided 
a total bandwidth of 4 GHz (after combining both IF bands). Total duration of the track 
(including integration time on source and calibrators) was ~6hr. IRAS 18276— 1431 was 
observed ±2.6 hr around its transit. 

The calibration of the data was performed using the MM A software package^. Data were 
gain calibrated in baseline-based mode using the quasar J1833— 210, which was observed at 
regular time intervals of ~ 20 minutes before and after our target. The quasar 3C273 and 
the planet Uranus were used as passband and flux calibrators, respectively. The flux of the 
gain calibrator J1833— 210 was estimated to be 2.2 Jy at 115 GHz. Flux cahbration errors 
could range up to 30%. 

Reconstruction of the maps from the visibilities was done using the Multichannel Image 
Reconstruction, Image Analysis and Display (MIRIAD) software. We Fourier transformed 
the measured visibilities with robust weighting, which is an optimized compromise between 
natural and uniform weighting, for ^^CO and natural weighting for ^^CO for S/N optimiza- 
tion. After that, data were cleaned and maps restored. The clean beam for our ^^CO (^^CO) 
maps has FWHM=10'.'9x7"l (FWHM=ll'.'5x8'.'7) and is oriented at PA=73° (PA=58°). 
The noise (la) in our ^^CO (^^CO) maps, as measured in channels with no signal, is 0.12 
(0.08) Jybeam"^. In our 3mm continuum maps cr=4mJybeam^^. The conversion factor 
from ^^CO (^^CO) surface brightness to temperature units is 1.19 (1.01) K per Jybeam"^. 
The coordinates of the tracking center in our maps are RA=18''30'"3K75, Dec=-14°28'56'.'28 
(J2000.0). 



3. Observational results 
3.1. Optical and NIR Imaging 

The NIR and optical images of IRAS 18276-1431 are presented in Figs. 1-4. IRAS 18276-1431 

displays a clear bipolar morphology, with two lobes oriented along PAf5il7°-20° separated 
by a dark waist. Our AO images have higher angular resolution, probe deeper through the 
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^MMA is written and maintained by the Caltech's Millimeter Interferometry Group. 
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dusty nebula, have larger dynamical range, and provide more structural details than the 
HST images. In particular, they show a sharp drop of intensity at the lobe edges, which 
are quite well defined. The lobes also have limb-brightened tips (or caps), suggesting that 
they are closed-end structures. These caps arc not visible in the optical images. The overall 
size and morphology of the two lobes are different, the SW lobe being shorter and narrower 
than the NE lobe. The morphological differences between the lobes are most prominent at 
the shortest wavelengths (optical and Xp-band images), where the SW lobe appears as a 
"two-finger" or A-hke structure in contrast with a more rounded NE lobe (Fig. 2 Right). 

The thickness of the dark equatorial waist between the lobes is larger in the optical 
images than in the NIR ones. This is because the optical depth of the material in the 
nebular equator obscuring the base of the lobes is expected to decrease with wavelength. 
Radial cuts of the NIR intensity along the lobes (Fig. 5) show a central minimum, where 
the extinction reaches a maximum value, and two adjacent absolute maxima (or peaks) 
separated by a distance 2/i=0'.'34 in the Lp and Mg images. (In the Xp-band image the 
separation between the peaks is slightly larger presumably due to a larger optical depth 
in Kp). The NIR optical depth along the line of sight at the peaks (i.e. at ±h from the 
central minimum) is expected to be <<1 given the sudden change in the slope of the radial 
profile (from a positive to a negative value): we start to see the base of the lobes through 
the waist. This probably indicates that the optical depth and, therefore, the column density 
of the waist decreases with latitude or, equivalently, with the distance above the equatorial 
plane. A density power law that decreases with the latitude is also inferred from modelling 
of polarimetric NIR data of IRAS 18276-1431 by Gledhill (2005). 

Our deep NIR images also reveal a faint (sensitivity limited) , elongated halo surrounding 
the lobes. The contribution of the extended PSF halo to this component is neghgible, less 
than ~l%-2%. The size of the halo is 4'.'5x3'.'4, as measured to a ~3(7 level over the sky 
background in our Kp image, and decreases with increasing wavelength, suggesting that 
the halo consists of scattered starlight. The decrease of size with wavelength may be also 
partially due to the lower sensitivity in the Lp- and Mg-band images. The major axis of 
the halo is oriented along PA~17°, similarly to the lobes as seen in the optical and NIR 
images. (In contrast, the extended PSF halo, which is observed around bright field stars, is 
elongated towards the guide star.) In the Lp image, there is faint emission extended in the 
direction roughly perpendicular to the lobes, along PA=105-110°. The outer radius of this 
sensitivity-limited structure is 0'.'6. 

Our Kp image has revealed two pairs of radial "searchlight beams" emerging from the 
nebula center and directed oppositely through the nebular center along PARi36° and PARi8° 
(Fig. 2). The angular width of the beams in the azimuthal direction is APAfslS". The 
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searchlight beams are not axially symmetric, the beams running along PA=36° being the 
brightest and the longest (observed up to ±2'.'8 from the nebula center). The north searchlight 
beam directed along PA=8° is much fainter than its southern counterpart. The south beams 
seem to be the continuation of the inner and brighter A-like feature observed in this lobe. 
There are several arc-like features superimposed on the searchlight beams (see also Fig. 6). 
Some of these arcs are also observed with a much smaller brightness in the regions of the 
halo between the beams. The separation between the arcs is in the range 0f!l8 to 0'/24 and 
some of them appear to intersect each other. The characteristic width of the arcs is 

The point where the beams cross at the nebula center (hereafter, beam center) docs 
not coincide with the center of the halo (Fig. 2). The latter has been determined by fitting 
elliptical isophotes^ to the halo in Kp band: we have only used the isophotes with interme- 
diate brightness, making sure that the main derived properties (center and shape) are not 
strongly affected by the brightness distribution of the lobes or bright field stars. The center 
of the beams has been estimated by taking azimuthal cuts with different center points along 
the polar axis. The center of the beams is defined to give the largest degree of symmetry 
in the width and position of the beams in the azimuthal cuts (we are assuming that the 
beams are intrinsically symmetric). The point where the beams cross at the nebula center 
and the center of the halo are separated by 53 mas along the direction of the lobes, which is 
larger than the expected errors in the location of the centers (errors < 1 pixel~9mas) . Partial 
illumination of the halo by the light scattered in the lobes (in addition to the starlight) could 
explain the offset between the two centers and the similar orientation of the halo to that of 
the lobes: the center of the halo is expected to move in the direction of the brightest lobe, 
i.e. the NE lobe. 



3.2. Photometry 

The optical and NIR photometry (and other observational details) are given in Table 
1. Total fiuxes have been estimated by integrating the surface brightness of all pixels within 
a polygonal area enclosing the emission from our source with signal 3a above the sky and 
leaving outside bright field objects. Our fluxes in the visible are signiflcantly smaller than 
those obtained from ground-based observations previously reported in the literature (see 
Fig. 7). This is because of the large apertures used in ground-based observations, which most 
likely included flux from the bright star FSl north-west of IRAS 18276— 1431 (Fig. 2). Our 
flux measurements in the NIR arc also smaller than the 2MASS fluxes but are in excellent 



•^Using the IRAF task stsdas. analysis. isophote. ellipse 
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agreement with the J and K photometry by Gledhill (2005) obtained using an eUiptical 
aperture including only flux from IRAS 18276—1431. The integrated flux in the NE lobe is 
larger than in the SE lobe by a factor 3-4 in the optical and 1.3-1.4 in the NIR (Table 1). 

3.3. Color maps 

We have obtained color images of IRAS 18276-1431 (F814W-F606W, L^-K^, M^-Kp, 
and Mg- Lp), which are shown in Figure 8. Registration of the HST optical and NIR images 
has been done using fleld stars to compute relative object shifts^. The estimated uncertainty 
in the image alignment is ~1 pixel (i.e., ~10 and ~23mas for the AO and HST images, 
respectively) . 

The reddest region in our F814W— F606W color map is not in the waist but in the base 
of the SW lobe, in a small region north of the relative maximum in the F814W and F606W 
surface brightness distribution. There is a color relative minimum at the waist, which is 
bluer than the base of both lobes. 

In our higher-resolution, better S/N ratio NIR color maps, the nebular structure is much 
more clearly delineated. The Lp—Kp and Mg—Kp color images show an hourglass structure 
with an overall point-like symmetry. The lobe walls and caps are quite well delineated in 
these images, in particular, in the Lp—Kp maps. The West (East) side of the North (South) 
lobe appears brighter. There is a relative color minimum (i.e. a relatively blue region) in the 
middle of the lobes, at a radial distance of 0'/3 from the center. As in the F814W— F606W 
color images, the Lp—Kp and Mg—Kp color images show a relative minimum in the waist. 
The elongation of the halo along PA=105-110° observed in the Lp-band images, is also 
noticeable in the Lp—Kp image with relatively blue colors. 

There are gradually increasing differences in the maps in going from Lp—Kp to M^—Kp 
and to Ms—Lp-. the lobe walls and caps smooth out and the waist reddens (relative to the 
lobes). In the Mg—Lp map the color reddens gradually from the lobe tips to the waist, 
where the nebula is reddest. These trends are mainly the result of the lower scattering and 
absorption optical depths at these long wavelengths. The waist in the Mg-Lp image shows a 
similar orientation to the more extended equatorial structure seen in the Lp-Kp, and Lp-band 
images. 



''Using the IRAF task images. immatch.imalign 
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3.4. Millimeter-wave observations 

Unresolved ^^CO J=l-0 emission is observed towards IRAS 18276—1431 within the LSR 
velocity range [+45:80] km s~^ (Figs. 9 and 10). The line profile is well fitted by a parabola 
centered at Vsys=62km s^^ (LSR) and total width (at zero intensity level) of 34 km s"^ The 
parabolic profile is characteristic of unresolved optically thick CO emission and yields an 
expansion velocity of the molecular envelope of IRAS 18276-1431 of T4xp=17km s^^ This 
value is comparable with, but slightly larger than, that derived from previous single-dish 
^^C0J=2-1 spectra and OH maser observations (Kxp=12-14km s~^; Heske et al., 1990; 
Bowers, Johnston & Spencer, 1983). The observed ^^CO profile has a prominent, narrow 
absorption feature in the range [+64:+69] km s~^(LSR), i.e., bluewards of Vgys, which is most 
likely due to interstellar absorption. For that reason, the integrated fiux over the line profile 
is expected to be a lower limit to the intrinsic ^^CO fiux from IRAS 18276— 1431. A better 
estimate of the intrinsic flux from the source is given by the integral of the parabolic flt to the 
line profile, which we have used for the determination of the molecular mass in the envelope 
of IRAS 18276-1431 (Table 2 and § 7.2). The ^^CO emission integrated over the fine profile 
is centered at RA=18^30"^3ffi67 Dec=-14°28'5r.'6 (J2000.0), with an error of ~0'.'8 both in 
RA and Dec. 

^^CO J=l-0 emission is detected above a 2a level in channels [+53:+58] km s~^ (LSR). 
Given its location, which coincides with the ^^CO and continuum source, the ^^CO fine most 
hkely arises in IRAS 18276— 1431 but the observed profile is affected by strong interstellar 
absorption at LSR velocities [+64:+69] km s~^ or even a broader velocity range, preventing 
us from obtaining an accurate characterization of the line profile. As for ^^CO, the measured 
^^CO line fiux is most likely a lower limit to the total fiux emitted by IRAS 18276— 1431. 
Parameters derived from our ^^CO and ^^CO maps are given in Table 2. 

An unresolved, continuum source is detected towards IRAS 18276— 1431 with a fiux of 
F2.6mm=15±4 mJy. The contamination of the continuum fiux by ^^CO 1-0 line emission is 
less than 5%. 



4. Extinction towards the central star 

The extinction along the line of sight towards the central star in IRAS 18276— 1431 
must be quite high not only in the visible but also at the observed infrared wavelengths 
since the illuminating star remains hidden even in the Mg band. We have estimated a lower 
limit to the optical depth (r) at 2.12, 3.76 and 4.67/xm from the ratio between the ob- 
served surface brightness (Jobs) and the expected unabsorbed intensity (/) per pixel from 
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the central star. Adopting values for the stellar effective temperature and luminosity of 
Teff=7000K and L=6AxW^Lq, a distance to IRAS 18276-1431 of i:)=3kpc (§ 1), and tak- 
ing into account the PSF in our images, we derive: /2.i2/^m=4.3xl0~^^, 13.76/^^=9.9x10"^'', 
and /4.e7^m=1.7x 10^^^ ergs^^ cm^^ per pixel. The non-detection of the star at our 
3cr level (Table 1) implies that the total (circumstellar plus interstellar) extinction towards 
the nebula center along the line of sight is T2.i2/im>12.6, r3.76^m>8.6, and T4.67/im>5.4. We 
have derived the interstellar extinction towards IRAS 18276— 1431, Ay=1.6±0.5mag, using 
the numerical algorithm provided by Hakkila et al. (1997). For a standard total-to-selective 
absorption ratio Ry=3.1, the above value of Ay yields the following estimates of the in- 
terstellar extinction: A2.2/im=0.3mag, A3.76^j„=0.16mag, and A4.67^ni=0.1. Adopting these 
values, the circumstellar optical depth towards the star in IRAS 18276—1431 is r2.i2>12.3, 
r3 76>8.4, and T4 67>5.3. These values are in agreement with the order of magnitude of the 
equatorial optical depth in the J band, rj^lO, inferred from modeling of the polarimetric 
data (Gledhill 2005), and the mean optical depth in the 10-25/Lim range, Tio-25^m^l, ob- 
tained from SED modeling (van der Veen et al., 1995; sec also §7.1), which indicates that 
the equatorial envelope around IRAS 18276—1431 is optically thick. Since the optical depth 
in the F-band (0.55/im) is typically a factor ~5 larger than in X-band (2.12//m) (for sihcate 
dust Ax oc A~^, see §7.1) the optical extinction towards the star must be ^y>60mag. This 
value is well above the lower limit to Ay derived directly from our HST F606W and F814W 
images following the procedure described above, To.606//m>15.4 and To.8i4//m>14.1. 



5. The illuminating source at the core of IRAS 18276—1431 

In this section we analyse the NIR colors of IRAS 18276—1431 and show that the lobes 
are illuminated partially by the starlight and partially by NIR light from a warm dust 
component, as has been also inferred for the bipolar PPN IRAS 16342-3814 (Sahai et al. 
2005). The presence of warm dust around the central star of IRAS 18276— 1431 is also 
deduced from the analysis of its full SED (Bedijn 1987 and §7.1). 

We analyse the observed NIR colors averaged over small apertures (of 0—Of!l) located 
in the middle of the lobes, where a relative minimum in Lp—Kp is observed. These regions 
(marked with circles in Fig. 8) are well-removed from the lobe edges and the dense central 
waist and, therefore, the optical depths to scattered light are expected to be smallest. In the 
NW- (SE-) lobe we find that the interstellar extinction corrected intensities (/^) for 2.1, 3.8 
and 4.7/im are in the ratio 1:3.8:7.8 (1:3.7:8.5) indicative of a very red source illuminating 
the lobes (note that the two lobes give very similar colors). The observed colors in these 
regions result from a combination of reddening due to internal extinction, and blueing due 
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to scattering, of light from the central illumination source. Although H2 emission has been 
detected in IRAS 18276—1431, line emission is not expected to be a significant contribu- 
tion to the emission in the NIR filters (Lowe et al., private communication — see also the 
featureless spectra of this object in the iJ, and L bands reported by Le Bertre et al., 
1989). Therefore, the observed intensity of the scattered continuum in the reflection lobes 
(corrected for interstellar extinction as in Sect. 4), /, is given by: 

/- /oe-"(l - e"^'-) (1) 

where, /q is the emission from the central source, r is the extinction (absorption+scattering) 
optical depth from the central source to the scattering region and from the scattering region 
to the observer (i.e. radial+tangential) , and rf^g is the hue of sight scattering optical depth 
at the observed region. The above expression is valid for scattering occurring predominantly 
in the lobe walls, which is a reasonable approximation for IRAS 18276—1431 (§6.1). 

We have modeled the K^-Lp and Mg—Lp colors (allowing ±20% errors) adopting the 
above formula and varying the temperature of the illuminating source. The absorption 
and scattering optical depths have been assumed to follow a power-law dependence with 
wavelength of the type ~A~^ and ~A~^, respectively, in the 2-5/im range, as expected for 
silicate dust (e.g. Draine & Lee 1984; Suh 1999). We find that the NIR colors measured 
cannot be reproduced if the only illuminating source is the central star, which is expected 
to have an effective temperature in the range Teff=4000- 10000 K. To explain the NIR colors, 
a warm (300-500 K) dust component is needed. In particular, the best fit of the NIR colors 
is obtained if the nebula is illuminated partially by the central star (for which we adopted 
refr=7000K) and partially by dust at 400 K. Our best fit model indicates that a significant 
fraction (~80%) of the central star's hght is absorbed by this warm dust cloud, which 
presumably covers a significant solid angle around the star. As shown in § 13, 400 K dust is 
expected to be located at 7-9xl0^^cm (45-60 AU) from the central star. Our best fit NIR 
color model also indicates that the light from both the star and the dust cloud is extinguished 
by T2.i2/im~3.7 in its way to, and through the lobes. 

6. Structural components in IRAS 18276-1431 

6.1. The bipolar lobes 

Our NIR color maps and images suggest that the lobes of IRAS 18276— 1431 are most 
likely cavities excavated in the AGB envelope (i.e., structures with dense walls and more 
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tenuous interiors), the light from the central source being reflected by dust in the lobe 
walls. This is also believed to be the typical situation in many post-AGB nebulae and, in 
particular for IRAS 18276— 1431, the OH maser spatio-kinematic distribution suggests the 
presence of bipolar outflows carving out the polar regions of the envelope (§ 1). Our Lp—Kp 
color maps indicate that the lobes are wider than they appear in the Kp image and display 
a certain degree of point-symmetry (see Fig. 11). Such a point-symmetry indicates that the 
east and west walls of the lobes have different densities, with the west (east) wall of the NE 
(SW) lobe being the densest, and/or are illuminated differently from the central star. The 
searchhght beam oriented at PA=36° seems to go through the least dense lobe wall. The 
orientation of the lobes on the plane of the sky as measured in the NIR color maps (PA~15°) 
is different than in the Kp image (PA~20°). The polarization angle measured by Gledhill 
(2005), PA=23°, is very similar to the average PA of the searchlight beams revealed by our 
AO images, and is probably a better representation of the direction along which the bipolar 
cavities were excavated. The different orientation of the lobes in the color maps is most likely 
the result of large optical depth through the west and east sides of the NE and SW lobes, 
respectively. Although the north and south lobes appear quite different in the Kp image, 
they are similar in the color maps, suggesting that the real structure of the two lobes is 
not that different and the A-like appearance of the SW lobe is mainly an illumination effect 
(presumably the same that produces the searchlight beams, which are similarily oriented). 

Scattering structures will in general tend to smooth out in images at long wavelengths 
because the scattering efficiency decreases as X~^. This trend is observed in our images. The 
limb-brightened appearance of the lobe tips or caps in our NIR images, which are best seen 
in the Lp-band image, suggests that the scattering optical depth through these regions is <1. 
The caps are not observed in the optical images, which is consistent with the larger optical 
depths through the caps at shorter wavelengths. The structured lobe caps smooth out in the 
Ms image most likely because the radial and tangential scattering optical depth through the 
caps is too low at 4.67/xm: most of the light from the central source that reaches the inner 
edge of the caps escapes in the polar direction and very little is scattered towards us. The 
radial cuts of the surface brightness along the lobes in the Kp and Lp images show that the 
caps move outward in going from Kp to Lp. The position of the lobe caps are very similar 
in the Mg and Lp images. This implies that the extinction in the polar direction through 
the caps at 2.12/im is larger than 1. In the Lp-band the lower extinction (<1) enables the 
photons to penetrate deeper the lobe caps. These values enable estimating a range for the 
optical depth through the lobe caps, 4<ry<7, adopting an extinction power-law ~ A~^. 
This yields a range of (2-4)xl0''cm~"^ for the mean density in the ~0'.'l-wide caps, assuming 
a gas column density to optical extinction ratio of N^/Ay ~3xl0^^cm~^ per magnitude, 
which is expected for silicate grains and a dust-to-gas mass ratio of 5=200 (§ 13) Taking into 



-14- 



account the size of the lobes and adopting a similar density throughout the lobe walls and 
caps, the total mass in the lobe walls must be ^10~^ Mq. The lower limit to the mass arises 
in the expected decrease of the density with the radial distance to the nebula center. 

The redder F814W-F606W colors of the SW lobe can be interpreted as an indication 
that the SW lobe is the farthest of the two, the larger reddening being the combined result 
of 1) the presence of more dust in the halo and the extended equatorial disk lying in front of 
the distant lobe, and 2) more efficient forward scattering (which would make the closer lobe 
appear brighter and bluer than the far one). The NE lobe is also the brightest in all the 
observed bands, which also favors its being the closest lobe. Note, however, that intrinsic 
differences in the dust density between the two lobes would affect their relative colors and 
brightness (e.g. Sanchez Contreras, Gil de Paz, & Sahai, 2004). 



6.2. The halo 

The halo of IRAS 18276—1431 most likely represents part of the circumstellar envelope 
ejected during the AGB phase that has not been strongly affected by the lobe shaping 
processes. The halo is probably spherical, the elongated shape being the result of the way 
the different regions of the halo are illuminated, in particular, of the lower optical depth 
along the poles and much larger extinction in the equatorial direction. The outer radius of 
the halo is >2'/8, as given by the length of the longest searchlight beams in the Xp-band. 
The radial profile of the halo surface brightness (Fig. 12) can give us information about the 
spatial distribution of the density in the AGB envelope and how the mass-loss rate that led 
to this component might have varied in the past. 

A cut of the Xp-band image along the nebular symmetry axis shows that beyond the 
tips of the lobes the halo radial profile varies with the radial distance to the nebula center 
like r^'^. A similar variation with r is observed in the Lp and Mg images. In the case of 
an optically thin scattering halo, the surface brightness varies proportionally to the light 
intensity dilution factor and to the line of sight scattering optical depth {t^^^) at each point. 
The light seen by the dust in the halo along the polar direction (PA=20°) emanates mainly 
from the central source, which is expected to be point-hke. Therefore, the light intensity 
dilution factor will vary as 1/r^. After correcting the observed flux distribution by this factor, 
we deduce that rfj^ oc in the halo. Since oc r^~" for a density power law oc r~", we 
conclude that a=3 in the halo. This result agrees with the dust density distribution ocr~^'^ 
deduced from SED modeling of IRAS 18276—1431 (van der Veen et al. 1995). The previous 
result is consistent with the AGB mass-loss rate increasing with time, assuming that the 
expansion velocity of the wind is time-invariant. 
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The NIR surface brightness radial profile along the equator shows a dependence 
for radii larger than 0'/6. This is consistent with the flux from these outer halo regions 
being mainly due to second (or higher) order scattering. In fact, these equatorial regions are 
mainly illuminated by the two bright scattering lobes rather than by direct starlight, given 
the very large optical depth along the equatorial direction [Ay > 60 mag). In this case, the 
light intensity dilution factor is expected to have a dependence with r shallower than that 
for a point-like illuminating source, probably of the type 1/r, which is the dilution of a flux 
field emanating from an infinitely long wire (see Alcolea al. 2000). After correcting by a 
~l/r dilution factor the observed flux distribution at radii >0'/6, we obtain riosocr"^ which 
is consistent with a density law pocr~^ as derived for the polar regions. 

The searchlight beams and arcs observed in the halo are most likely illuminated parts 
of the AGB shell viewed in scattered hght. Discussion on the very uncertain origin of these 
features is deferred to Section 8. 



6.3. The waist 

The two lobes of IRAS 18276— 1431 are separated by an equatorial dust waist/torus 
behind which the central star is hidden. OH maser emission is observed to arise in this 
equatorial torus over a range of angular radii of 0'/25 and 0'.'85 (Bains et al., 2003, § 1). The 
spatial coincidence between the NIR dark waist and the OH maser emitting region is expected 
since the OH masers are pumped by FIR radiation from dust and are more likely to occur 
where the density is highest. The faint, equatorial structure seen in the Lp-band image with 
outer radius Rout ~0'.'6 (§3.1) probably corresponds to outer parts of the equatorial torus, 
which are seen through multiply scattered light: a fraction of the photons scattered in the 
lobes are directed towards the equatorial plane and from there scattered again toward us. 
Although equatorial disks/tori are usually inferred from the presence of dark lanes separating 
the bipolar lobes, there are other PPNs were these equatorial structures are observed in 
scattered hght, e.g. IRAS 04296+3429 (Sahai 1999), IRAS 17106-3046 (Kwok et al. 2000), 
and IRAS 17347-3139 (Sanchez Contreras et al. 2006). The relatively blue colors of the 
waist in the L^,—Kp and Ms—Kp color maps (§3.3) are consistent with multiple scattering 
illuminating with blue photons these equatorial regions, which otherwise would be expected 
to show the reddest colors. Since scattering is much less efficient at long wavelengths, the 
inner parts of the waist are best delineated in the Mg—Lp color map, where the waist displays 
the reddest colors, revealing the regions with the largest extinction along the line of sight. 

As shown in Section 4, the optical depth of the torus towards the central star along the 
line of sight is 74y>60mag. This value yields a column density of gas in front of the star of 
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Nd> 1.8x10^^ cm~^, for a gas column density to optical extinction ratio of Nd/Ay ~3xl0^^ cm~^ 
per magnitude (see §6.1). Adopting an outer radius for the torus of Rout=0'.'85, as given by 
OH maser emission mapping (see above), and assuming that the inner radius of the torus is 
much smaller than Rout (as indicated by SED modeling, §7.1), the mean H2 density in the 
torus towards the central star along the line of sight is >5xlO^[D/3kpc]~^ cm~^. 

An order of magnitude for the total mass in the waist enclosed within a radius of 
-R<0'.'85 and an height of h<0'!l7^ can be roughly estimated adopting a given density law 
as a function of the height above the equatorial plane, z. In the following, we consider two 
extreme cases for p{z). For a constant density throughout the torus, /)(2;)=constant, the total 
mass in this component would be >4 [i^/3kpc]^ Mq. A more realistic density law is one in 
which the density decreases with increasing z (§3.1). For example, for a steep exponential 
density law of the type p{z)—p{0) exp{—zlnx/h) where x=p(0)/p(/i), the mass in the torus 
would be >4x(l-l/x)/lnx [D/3kpc]^M0. Our NIR images suggest that the optical depth 
through the waist at h is probably <1 in the 2-4/xm wavelength range, since we start to 
see through the base of the lobes (§3.1). Since the extinction in the optical is typically a 
factor 5 larger than in the i^p-band, we estimate that Av{h)<5. Therefore, a reasonable 
lower limit to x (ocNd(0)/A'd(/i)ocAy(0)/Av'(/i)) may be ^10, which yields a total mass of 
~ 1 [£'/3kpc]^ Mq. The mass derived is very uncertain since the total extent of, and density 
distribution in, the torus as well as the geometry of the innermost nebular regions, which 
determines the way the nebula is illuminated, are unknown. Therefore the above figure has 
to be taken only as an indication that the mass of this equatorial component could be quite 
large, in particular, it could represent a significant fraction of the total mass in the envelope. 
This is also observed in other post-AGB nebulae, for example, in OH 231.8+4.2, which has 
a dense (~3xl0^cm~^) equatorial waist that contains ~70% of the total envelope mass as 
derived from high-angular resolution CO mapping (Bujarrabal et al., 2002; Alcolea et al., 
2001). 

For a silicate+graphitc mixture in the envelope of IRAS 18276—1431, which cannot be 
completely ruled out (sec §13), the density and the mass derived for the torus would be 
smaller. For example, the standard ISM graphite and silicate mixture with a dust-to-gas 
mass ratio 5=200 yields N^/Ay ~1.5xl0^^ cm~^ per magnitude. For this lower value of the 
Nd/Ay ratio, the mean density in the torus would be >2xl0^ [D/3kpc]~^ cm~^, leading to 
a mass of ~(0.2-0.07) [D/3kpc]^M0 for the two extreme (constant and exponential) density 
laws considered above. 



^Projection effects can be ignored since the lobes are very close to the plane of the sky and, therefore, 
the torus is expected to be nearly edge-on (§ 1). 
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7. Mass, mass-loss, and dust temperature 

7.1. SED modeling 

We have derived the physical parameters of the dust envelope of IRAS 18276—1431 by 
using a radiative transfer code to fit its SED at wavelengths A>5/xm, the region of the 
spectrum where dust thermal emission is the main contributor to the observed flux (Bedijn 
1987). At A<2/xm the dominant contribution to the flux is the scattered light from the 
central source. In the range 2 to 5/im, a combination of these two emission mechanisms is 
expected. Calculations were performed with the photoionization code CLOUDY (version 
05.07) last described by Ferland et al. (1998). CLOUDY has a dust code built in that 
solves the time-independent radiative transfer equation under the constraint of radiative 
equilibrium. Details of the grain model are given in van Hoof et al. (2004). The code 
makes a comprehensive, detailed treatment of the grain physics that includes, e.g., quantum 
and stochastic heating, heating of the grains by Lya photons, and the photo-electric effect. 
CLOUDY does not work under the optically thin approximation, but calculates the thermal 
dust continuum for optically thick envelopes. 

We have assumed that the envelope is spherically symmetric. Despite the simplicity of 
this assumption, the code is expected to provide reasonable constraints to the mean proper- 
ties of the dust envelope (mass and temperature). As shown by e.g. Ueta & Meixner (2003), 
asymmetries in the density distribution would mostly affect the optical part of the SED but 
not the MlR/FlR, which is the spectral region we are modeling. Scattering is assumed to be 
negligible in our model. As discussed in detail by Bedijn (1987), this is a reasonable approx- 
imation at A>2/i,m for dust shells with intermediate optical depths (i.e. r9,7^m~4), which 
is the case of IRAS 18276— 1431. Even for very optically thick envelopes, with r9.7^^~30, 
inclusion of scattering does not change significantly the SED in the wavelength range that 
we are modeling. 

In order to calculate the SED of IRAS 18276—1431 we have adopted a series of input 
parameters. We have assumed that the central star has a blackbody spectrum with an 
effective temperature of Teff=7000K and bolometric luminosity of L=64OOL0(§ 1). We have 
adopted astronomical silicate grains as suggested by the presence of a weak 9.7/im silicate 
feature in the IR spectrum of IRAS 18276-1431 (Le Bertre et al. 1989; Bedijn 1987). The 
optical constants for astronomical silicates are from the calculations of Martin and Rouleau 
(1990), which extend the work of Draine and Lee (1984) to ionizing energies where the grains 
are strongly absorbing. The high IR polarization measured in IRAS 18276—1431 suggests a 
population of small dust grains, with radii a<0.3/im (Gledhill 2005). Therefore, we have used 
the grain size distribution adopted by these authors, that is, with minimum and maximum 
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grain radii of ainin=0.01//m and ainax=0.3//m and a power law size index of —5.5. CLOUDY 
uses a spherical Mie code to generate the grain opacity as a function of wavelength for a 
given size distribution and grain material (van Hoof et al., 2004). 

The other input parameters needed by the code to calculate the SED are the inner and 
outer radius of the envelope, and the density distribution. For the density, we have adopted 
a radial power-law p oc r~", and a dust-to-gas mass ratio constant throughout the envelope 
with a value of 5=200, which is typical of post-AGB 0-rich envelopes. As already suggested 
by Bedijn (1987), a two-shell structure with warm and cool dust is needed to reproduce 
the SED in IRAS 18276— 1431. The inner, warm dust shell would represent a recent low 
mass-loss rate wind that started after the cessation of the heavy mass-loss at the end of the 
AGB. The former period of heavy mass-loss led to the outer, cooler envelope. In the outer 
shell, we have adopted a—?>, as indicated by the radial profile of the halo in our NIR images. 
For the inner shell we have used a— 2, which is expected for a time-invariant mass- loss rate 
at a constant expansion velocity. The free parameters of oTir model arc: the inner radius of 
the inner shell, -Rin,i, the inner radius of the outer shell, i?in,2 (which is also the outer radius 
of the inner shell), the outer radius of the outer shell, -Rout, and the density (or equivalently 
mass-loss rate for Vcxp=17km s~^) at the inner radius of each shell. These parameters have 
been varied until a good match of the observations is obtained. The best-fit parameters have 
been determined calculating and minimizing a function: 

= (i^A,obs - F^Y/al (2) 

where a\ is the error of the observed flux Fx,obs, which has been assumed to be of 20% 
for all wavelengths in the range modeled (i.e. 5-100//m). 

The parameters describing the dust envelope for the best-fit model are given in Table 3. 
The resulting synthetic SED together with the radial dust temperature distribution are 
plotted in Fig 13. In this figure we also show the gas density distribution for the best-fit 
model for a constant gas-to-dust mass ratio throughout the envelope (see above) and the 
kinetic temperature of the molecular gas predicted by CLOUDY after solving the full physics 
(heating and cooling processes) for the derived envelope parameters. 

The total mass in the envelope of IRAS 18276—1431 enclosed within an outer radius 
of i?out=l-25xl0^^cm is ~O.95M0, most of it is in the outer, cooler shell. Larger values of 
-Rout reproduce equally well the SED of IRAS 18276—1431, therefore, the value 0.95Mq may 
be a lower limit to the total mass. (For an outer radius of lO^^cm the total mass of the 
envelope would be 1.9M0.) Although the best-fit SED model is obtained for an inner radius 
of i?in,i=6.3xl0^'^cm, values in the range i?in,i=[5-8] xlO^^cm lead to acceptable results. The 
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optical depth is Ay~4mag, which is in good agreement with the optical depth of the 9.7//m 
sihcate feature derived by Bedijn (1987): note that the opacity of sihcate grains as a function 
of wavelength is such that Ay^Tg^T^m- The extinction derived from SED modeling is much 
smaller than the extinction towards the central star along the line of sight (>60 mag) derived 
from our NIR images. This is expected since SED modeling is done assuming spherical 
symmetry and, therefore, only provides average values for the envelope parameters. The 
relatively low value of Ay derived from SED modeling compared with Av>60meig indicates 
a large equator-to-pole density constrast, consistent with the optical and NIR images. The 
dust temperature in the inner shell can be approximated by a radial power-law of the type 
7dust('") ^ r~°'^, except for the innermost regions where the temperature increases with a 
steeper slope. In the outer shell, the dust temperature decreases outwards following a slightly 
shallower radial power-law, T(iust('") r"*^-^^. Average gas densities in the inner and outer 
shell are in the range f«10''-10 cm" -3 and cm , respectively. Densities could be a 

factor >15 larger in the equator and significantly smaller along the polar direction. 

Our model clearly underestimates the observed fluxes in the submillimeter and millime- 
ter wavelength ranges. Adopting a larger value for i?out has little effect on the flux at long 
wavelengths given the steep density law in the outer shell, therefore, a larger i?out cannot 
explain the emission excess. The emission excess could be partially accounted for by a dust 
emissivity function shallower than that assumed in our model, Qx oc A~^, at wavelengths 
>100/xm. SED modeling of a number of AGB and post-AGB objects in the submm and mm 
range suggest that the absorption efficiency of the dust in 0-rich envelopes follows indeed a 
A"'^-^ to A"^ law (Knapp et al. 1993 ; Knapp et al. 1994). However, from theoretical studies 
and laboratory works, an exponent <1.25 is only justified with difficulties (Giirtler et al.l996, 
and references therein). The submm and mm fiux excess can also be explained if there is 
an additional component of larger grains, as it has been inferred for a number of post-AGB 
objects (e.g. Walmslcy et al. 1991, Jura et al. 1987, Sanchez Contreras et al. 1998, Sahai 
et al. 2006). The SED of IRAS 18276-1431 in the 800/im-2.6mm range is very well fitted 
by that of a black-body in the Rayleigh- Jeans limit, i.e. fi, oc i/^ (Fig. 13). Since grains 
are expected to radiate like blackbodies at wavelengths A<27ra, an additional population 
of grains with radii a>400//m would explain the observed spectral index of the continuum 
at long wavelengths. The temperature of this large-grain-component must be in the range 
~150-20K: if Tdust>150K then we would obtain too much flux at A<100/im; also Td must 
be >20K since the spectral index in the 800;um-2.6mm is consistent with black-body emis- 
sion in the R-J limit. Since the dust model built in CLOUDY does not allow in its current 
version grains larger than 5/xm, we have estimated the mass of large grains needed to ac- 
count for the mm continnuum flux as explained in Sanchez Contreras et al. (1998). We have 
adopted a mean grain density and radius of p—3 gcm~^ and a > 400/xm, respectively. As- 
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suming that the dust emissivity at A<2.6mm is Q\^l, we estimate that a total mass of dust 
of Afdust>1.6xl0-4 [D/3kpc]2MQ, for Tdust=150K, and Md„st> 1.6x10-3 [D/Skpc^MQ, for 
Tdust=20K, is contained in this big grain dust population. Adopting a gas-to-dust mass 
ratio of (5=200®, the total gas mass associated with this large-grain-component is M~(0.03- 
0.3) [D/3 kpc]^ [5/2OO]M0. The location of this large-grain-component is unknown, since 
the mm-continuum emission is unresolved in our ~8"-resolution maps. Cold (20-30 K), big 
grains have been inferred to be located in the dusty equatorial regions of some PPNs, where 
the grains are shielded against the starhght due to the large extinction along the equator 
(e.g. IRAS 22036-1-5306, Sahai et al., 2006). The temperature of the large-dust-component 
in IRAS 18276—1431 is not well constraint and, therefore, there is not a sufficiently strong 
reason to claim that the grains are in the dark (cold) waist. Big grains are also believed 
to grow in long-lived, orbiting disks (e.g. Red Rectangle, Jura et al. 1997). High-angular 
resolution mapping of the sub-mm or mm continuum emission and/or measurements of the 
continuum in the A~300-500yLim range are needed to unveil the location of the big grains in 
IRAS 18276-1431. 

Finally, our model predicts a 9.7/im-silicate feature that is slightly stronger than ob- 
served, however, we have not found any models that simultaneously produce a weaker silicate 
feature and provide a good fit of the SED. This may suggest that there is a certain mixture 
of grains with different composition in IRAS 18276—1431. 

7.2. The mass from CO data 

As for most PPNs, the CO emission most likely arises in an expanding envelope that 
is the remnant of the ACB envelope and that is probably the counterpart to the NIR halo 
and waist. We have estimated the mass of this component from the ^^CO and ^^CO 1-0 flux. 
We have assumed that both hnes are thermahzed, i.e. the CO excitation is described by a 
rotational temperature that it is similar to the kinetic temperature, Tkin- The thermalization 
assumption is a reasonable approximation for most PPNs and, in particular, is expected 
to be vahd for IRAS 18276— 1431, which has nebular densities larger than the COJ=1-0 
critical density (wlO^ cm^^) in most regions of the CO-emitting envelope probed by our maps 
(i?co<5") as shown in Fig. 13. We have used the mean kinetic temperature of H2 predicted 
by CLOUDY for our best-fit model (Fig. 13) to estimate the molecular mass from CO. We 
have adopted a mean temperature of Tkin=110K and a relative ^^CO-to-H2 abundance of 



^In the absence of reliable estimates of S for the large-grain component in circumstellar envelopes, we 
adopt the same value as for the 'normal' dust component. 
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2xl0~^ and ^^CO-to-H2 abundance of 10~^, which are common values in PPNs (Bujarrabal 

et al. 1997, 2001, and references therein), and the well known molecular parameters of 
CO (partition function, Einstein coefficients, etc; see e.g. Rohlfs and Wilson, 2000). We 
derive a total mass of molecular gas in IRAS 18276— 1431 of >0.28Mq from ^^CO (using 
the integrated flux of the parabolic fit to the line profile) and >O.38M0 from ^'^CO (Table 
2). The lower limit to the mass derived from ^^CO is because this line is optically thick 
as indicated by its parabolic profile. Additional confirmation is given by the intensity ratio 
between the ^^CO J=2-l line, which was observed by Heske et al. (1990) with the IRAM 
30 m radiotelescope (beam=21") with a main beam temperature of Tmh—^ K, and the J=l-0 
line, for which we measure Tmb^l-SK with a 10'.'9x7'/l beam. This yields a beam-corrected 
J=2— 1-to- J=l— line ratio of ~1, which indicates that both lines are optically thick through 
most regions of the CO-emitting envelope. The lower limit to the mass derived from the 
^^CO line arises because of the strong interstellar absorption of the circumstellar emission 
and maybe because of a large optical depth in this transition as well. The fact that the 
mass derived from ^^CO is smaller than that derived from ^^CO, which is expected to be less 
opaque, is in agreement with the -"^^CO J=l-0 transition being optically thick. 

The lower limit to the mass derived from CO is in agreement with the total mass 
obtained from SED modeling, ^IMq. The total mass enclosed within the upper limit to 
the outer radius of the CO-emitting nebula probed by our maps, i.e. within 5" (2xl0^^cm), 
is 1.2 M0, which is also consistent with the lower limit to the mass obtained from CO. 

We note that values of Tkin of the order of ^10^ K arc expected and observed in young 
PPNs with dense envelopes (e.g. CRL 618, Sanchez Contreras et al. 2004, see also Hrivnak & 
Bieging 2005). Values of Tkin as lows as 15-30 K are also observed, however, in a number of 
PPNs (Bujarrabal et al., 2001). For a mean value of the kinetic temperature of 50 K, which 
is the average of Tkin over volume for IRAS 18276—1431 (Fig. 13), the lower hmit to the mass 
derived from ^^CO is ~0.2Mp,. 



8. Formation and evolution of IRAS 18276-1431 

The analysis of our high-angular resolution NIR and optical images, and the SED of 
IRAS 18276—1431 enable us to trace back in time the evolution of this object and to attempt 
understanding how the different structural components revealed by our data were formed. 
From the values for the outer and inner radius of the outer dust shell of IRAS 18276—1431 
(-Rout^l -25x10^^ cm and i?in,2=l-6xl0^^ cm, Table3) and the expansion velocity derived 
from CO, \4xp=17km s~^, we infer that the heavy mass-loss process that led to this compo- 
nent started at least ~2300yr ago and ceased ~300yr ago. During that time the mass-loss 
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rate increased from M~2xlO~^M0 yr"^ to 2xlO~^M0 yr"^. This heavy, time-increasing 
mass-loss-rate wind, which is expected at the end of the AGB phase, was succeeded in 
IRAS 18276— 1431 by a lower mass-loss rate wind (M^SxlO'^M© yr^^), which could have 
stopped, or been replaced with an even more tenuous wind recently {^10-15 jt ago). The 
presence of bipolar lobes in IRAS 18276—1431 may imply large expansion velocities in some 
regions of the envelope (mainly along the polar axis) . Although most of the dust is expected 
to be located in the massive AGB envelope, some dust fraction is located in shock-accelerated 
regions (e.g., the lobe walls), which are expected to have a complex and unknown kinematics. 
Therefore, the ages derived above assuming T4xp=17km s~^ are uncertain, especially for the 
inner dust shell. 

As proposed for most PPNs (§1), the lobes of IRAS 18276— 1431 have been probably 
excavated by interaction of fast, bipolar (jet-like) winds with the AGB envelope. The jet- 
launching mechanisms which are believed to operate in the early post-AGB or late AGB 
phase are still the subject of a vigorous debate. Prom the deprojected total length of the 
lobes (0'.'5~2.2xl0^^ cm) and adopting an expansion velocity of ftilOOkm s~^, which is typical 
of bipolar post-AGB outflows (e.g. Bujarrabal et al. 2001), we estimate a value of ~65yr 
for the kinematical age of the lobes. As shown by our NIR color maps, the lobes display a 
certain degree of point-symmetry. Such a symmetry is found in many post-AGB objects and 
is often interpreted as a signature of precessing jet-like winds (e.g. Sahai et al. 2005). 

The origin of the equatorial waist / torus is very uncertain. Given the large extent of the 
torus compared to the bipolar lobes, we believe that the waist is unlikely the result from 
the interaction of fast jets with an initially isotropic AGB wind. If the waist is just what 
is left of the AGB CSE after the jet-l-'AGB envelope' interaction, then we should observe 
a spherical envelope with a central biconical cavity that is as large as the optical lobes. In 
contrast, the dense material in the waist obscuring the central star and the base of the lobes 
is constrained to the equatorial region forming a torus that extends up to a radius of >0'.'6 
(>2.7xl0^^cm), as derived from our Lp-band images, and maybe of >0'.'85 (>3.7x 10^^ cm), 
as derived from OH maser emission mapping (§ 1). It may be possible for shocks developed 
in the jet+'AGB envelope' interaction to move sideways and backwards compressing the gas 
in the equatorial plane and shaping the torus. However, since these shocks are expected to 
be much weaker and slower than the forward shocks which produce the lobes, they can form, 
at best, a torus whose radius is much smaller than the length of the lobes, contrary to what 
we observe. Alternatively, the waist of IRAS 18276—1431 could be the result of a period 
of enhanced mass-loss concentrated in the equatorial plane. However, what induces such a 
heavy and asymmetrical mass-loss is unknown, especially considering that the mass in the 
disk, although very poorly determined, could be a signiflcant fraction of the total mass in 
the envelope (§6.3). Adopting an expansion velocity of Vexp=17km s~^, we estimate that 
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the kinematical age of the torus is >700yr. 

Searchlight beams and arcs hke those observed in IRAS 18276—1431 are two of the 
most fascinating, yet poorly understood features revealed by recent high-angular resolution, 
high sensitivity images of AGB and post-AGB objects. Beams and arcs are most likely 
illuminated portions of the envelope viewed in scattered light. As proposed by Sahai et 
al. (1998), searchlight beams may result from starlight escaping through holes/cavities in 
an inner dust cocoon. The agent excavating the holes is unknown but it may be the same 
collimated (precessing) jets responsible for carving out the bipolar lobes. Morris (1981) 
had originally proposed another, more simple, scenario to explain the structure of bipolar 
reflection nebulae including searchlight beams, which were observed as small protrusion 
on the reflection lobes of some PPNs in low-angular resolution images. In this scenario, 
which was also adopted in the scattering models by, e.g., Yusef-Zadeh, Morris, & White 
(1984) and Latter et al. (1993), these protrusions result from a density law that declines 
with latitude above the equatorial plane if, at some latitude, the radial optical depth to 
starlight falls rapidly below unity. As extensively discussed by Sahai et al. (1998), the 
latitude-dependent density model, however, was not able to explain satisfactorily the detailed 
morphology of the searchlight beams in the PPN CRL 2688 as observed with high-angular 
resolution. Some of the objections to such model argued by Sahai et al. (1998) also apply 
to IRAS 18276—1431. In particular, our images do not show any evidence of systematic 
ellipticity in the arcs as would be expected in Morris' scenario. Moreover, the latitude- 
dependent density scenario does not account for the sharp edges and relatively large width 
of the searchlight beams in IRAS 18276-1431 (APARil5°; §3.1). Nor does it account for 
the point-symmetry of the searchlight beams as observed in om T^^p-band image (§3.1), 
which would require a point-symmetric density distribution at a spatial scale as large as the 
extended halo, where the beams are observed up to a radial distance of ±2'.'8. Such a density 
law would be very difflculy to support since it would imply a point-symmetric mass loss 
during the AGB phase. As far as we know, searchlight beams are also observed in the PPNs: 
CRL2688 (Sahai et al.l998), IRAS 20028+3910 (Hrivnak et al. 2001), IRAS 17150-3224 
(Su et al.2003), and IRAS 17245-3951 (Hrivnak et al. 1999), raising to flve the number 
of PPNs with searchhght beams. The increasing number of objects with beams suggest 
that these are not extremely rare features, and, therefore, the dust distribution required to 
explain the beams must be relatively common at least during some time in the evolution of 
these objects. It is worth noting that the 5 objects listed above display bipolar morphology, 
have dense, equatorial waists, and are seen nearly edge-on (IRAS 20028+3910 appears to 
be at an intermediate orientation, Su et al. 2001). The presence of a nearly edge-on, dark 
waist may then be a necessary condition for the observation of the beams. This is in part 
expected because scattering features will be necessarily more difficult to observe in objects 
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with very bright cores outshining the nebula. Another interesting aspect is that beams are 
not observed in PNs or AGB envelopes, but only in young PPNs. The lack of beams in PNs 
may be an observational effect, since their bright central stars are usually directly visible 
and the scattering nebula is very diluted. The lack of beams in AGB envelopes is consistent 
with the idea that the winds clearing out the holes in the cocoon start after the star has 
left the AGB. The spectral types of the searchlight beam's PPNs are concentrated in the 
F-G range. This may be partially due to destruction of the dust cocoon by the high- velocity 
outflow responsible for carving out the holes and/or by the increasing stellar UV field as 
the star evolves towards hotter spectral types. Finally, the presence/absence of searchlight 
beams does not seem to be strongly correlated with the chemistry, since beams are observed 
both in C- and 0-rich objects, however a larger sample of searchlight beam PPNs is needed 
to reach reliable conclusions. 

Arcs and rings are also getting more and more common in PPNs. Up to date arcs have 
been detected in 7 PPNs (including IRAS 18276-1431), 4 PNs, and 1 AGB envelope (see 
summary of arcs/rings in post- AGB objects by Su 2004). The binary-induced spiral outflow 
model originally proposed by Mastrodemo and Morris (1999) produces arc-like features very 
similar to those observed in IRAS 18276—1431 for a system viewed from near the equatorial 
plane and taking into account illumination effects (such as very large extinction along the 
equator by the dusty waist). In this model (see Fig. 14), the arcs disappear along the 
polar axis and alternate in radius on each side of the bipolar axis. The alternating arcs 
spaced along the polar axis are best seen if the angular momentum axis of the system is 
perpendicular to the line of sight. From our point of view, the most attractive aspect of the 
binary- induced spiral outflow model is the ability to explain the non-concentric, intersecting 
arcs of IRAS 18276— 1431. This model has also been proved to explain the spiral pattern 
recently discovered in the extended envelope of the PN AFGL 3068 (Mauron & Huggins, 
2006; Morris et al., in preparation). A possible argument against this model is that predicts 
arcs separated at quite regular intervals, set by the orbital period, whereas the time spacing 
of the arcs in IRAS 18276-1431 ranges between 190 and 250yr for V;xp=17km s^^ Models 
involving magnetic activity cycles (e.g. Soker 2000, Garcia-Segura et al. 2001) appear as an 
interesting alternative possibility to explain non-regularly spaced arcs. However, it is not 
known whether or not magnetic periodic variations can happen on a scale of hundreds of 
years, as required to explain the arc spacing in most PPNs (Su 2004). 
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9. Summciry 

We present and analyse NGS-AO images in the JCp, Lp, and Mg filters of the PPN 
IRAS 18276—1431 obtained with the Keckll 10m telescope. Supporting optical F606W and 
F814W HST images as well as interferometric observations of the ^^CO (J=l-0), ^^CO ( J=l- 
0), and 2.6 mm continuum emission with OVRO are also presented. Our main findings are 
summarized as follows: 

- Our AO images reveal the structure of the dust envelope around IRAS 18276—1431 
with O'.'06-O'.'l-resolution. The envelope of IRAS 18276—1431 displays a clear bipolar 
morphology with two lobes separated by a dark waist and surrounded by a faint, 
sensitivity-limited halo. Our i^Tp-band image reveals two pairs of radial "searchlight 
beams" emerging from the nebula center and several intersecting, arc-like features. 

- The lobes are probably cavities with dense walls and caps excavated in the AGB 
envelope presumably by tenuous bipolar winds. The point-symmetry displayed by the 
nebula may suggest that these bipolar winds are precessing. Although quite uncertain, 
our estimate of the total mass in the lobe walls is >10"^Mq. For an adopted lobe 
expansion velocity of 100 km s^^, which is typical of post- AGB outfiows, we estimate 
a kinematical age ~65yr for this component. 

- The halo most likely represents the remnant of the AGB CSE that has not been 
strongly affected by the process that produces and shapes the lobes. The density in 
the NIR halo follows a radial power-law ocr~^, which is consistent with a mass-loss 
rate increasing with time during the latest stages of the AGB. The outer radius of the 
halo is >1.25xl0^'''cm, which implies that the halo was ejected >2300yr ago, for an 
expansion velocity of yexp=17km (see below). 

- The illuminating star remains hidden behind the dark waist even in the Mg-band 
images. Such a waist most likely represents a dense equatorial torus, also traced by 
OH maser emission, that could have resulted from equatorially enhanced mass loss. The 
circumstellar optical depth towards the star along the line of sight (i.e., through the 
torus) derived from our NIR images is r2.i2>12.3, T3.76>8.4, and T4.67>5.3, which implies 
a large extinction in the optical, Ay>60mag, and a mean density >5xl0''' cm~^. The 
mass in the equatorial torus is very uncertain, but it could represent a large fraction 
of the total mass of the envelope. 

- The searchlight beams may result from a preferred illumination of the halo through 
"holes" excavated in an inner dust cocoon, maybe by the same agent responsible for 
carving out the bipolar lobes. 
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- The intersecting arcs, which are spaced at time intervals ranging between 190 and 
250 yr, are reasonably well explained by the binary-induced spiral-outflow model orig- 
inally described by Mastrodemos & Morris (1999), which produces arc-like features 
very similar to those observed in IRAS 18276—1431 for a system viewed nearly edge-on 
and taking into account illumination effects (Morris et al, in preparation). 

- The analysis of the NIR colors indicates that the scattering nebula is partially illumi- 
nated by the central star (with Teff =4000-10,000 K) and partially by a warm (~300- 
500 K) dust component. The presence of warm dust is confirmed by modeling of its 
full SED. 

- We have modeled the 5-100/im wavelength range of IRAS 18276—1431 using the dust 
radiative transfer model built in CLOUDY. Based on the detection of a weak 9.7/im- 
silicate feature, astronomical silcates have been used in our model. The grain size distri- 
bution used is the same as that adopted by Gledhill (2005) to explain NIR polarimetry 
data. The central star has been assumed to emit like a Tejf=7000K, L=6400Lq black 
body at a distance D=3 kpc. 

- SED modeling suggests a two-shell structure: 1) an outer dust shell with inner and 
outer radius i?in~1.6x lO^^cm and i?out^l-25x 10^''cm, dust temperature Td~105-50K, 
and mass-loss rate that increased from M~2xlO~^M0 yr^-*^ to 2xlO~^M0 yr~^ in 
~2000yr (adopting an expansion velocity VI,xp=17km s~^, see below); and 2) an inner 
sheU with i?i„~6.3xl0i4 cm, Tdust~500-105 K, and M~3x10-5Mq yr'^. 

- The mass of the envelope enclosed within i?out=l-25xl0^^cm derived from SED mod- 
eling is ^IMq, for a dust-to-gas mass ratio of 5=200. 

- The molecular envelope of IRAS 18276—1431 is unresolved in our ~8"-resolution ^^CO 
maps. Our CO data indicate that the bulk of the CO-emitting nebula is expanding at 
low velocity, T4xp=17km s~^. Adopting a mean kinetic temperature in the envelope of 
7kin=110K, as predicted by CLOUDY for the best-fit model to its SED, we derive a 
mass of M>O.38M0 for the molecular envelope. 

- Our best-fit SED model underestimates the observed sub-mm and mm fiux. This excess 
can be accounted for by an additional population of big grains (radius a>0.4mm) with 
temperature and mass in the range Tdust=150-20K and Mdust=(0.16-1.6) xlO~^M0, 
respectively. The location of this large-grain-component is unknown, since the 2.6mm 
continuum emission is unresolved in our maps. One possibility is that big grains are 
located in the dense equatorial waist. 
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Tabic 1: AO Sz HST observations Sz pliotoiiictry of IRAS 18276 — 1 131 



filter 


Ao 


AA 


time 


FWHM 


rms (cr) 


-^max 


Flux 


NE-to-SW 




(//m) 


(//m) 


(sec) 


(mas) 


(mJy/pix) 


(mJy/pix) 


(mJy) 


Flux ratio 


F606W^ 


0.61 


0.123 


200 


138 


4x10"^ 


3.8x10-3 


0.9 


4 


F814W^ 


0.81 


0.149 


120 


138 


1x10-5 


1.1x10-2 


3.0 


3 


K ^ 


2.12 


0.351 


270 


60-65^= 


8x10-6 


4.3x10-2 


96 


1.3 


T b 

Lip 


3.78 


0.700 


225 


80-85^= 


3x10-^ 


0.29 


385 


1.4 




4.67 


0.241 


180 


95-100^ 


2x10-3 


0.70 


781 


1.4 



°HST subsamplcd images: FoV=35"x35"; 22.8 mas pixer^ 
''Keck AO images: FoV=10"xlO"; 9.94 mas pixel 

■^Range determined from PSF fitting to several field stars at the same radial distance from the reference star 
than our target. 



Table 2: CO line parameters and derived molecular mass 





i^CO 1-0 


1300 1-0 


Fsys (km s-^) 


62 


a 


Vcxp (km s-^) 


17 


a 


Ico (Jybeam-^) 


1.5±0.3 


0.32±0.08 


/ Ico dV dn (Jykm s-^) 


27±3 (34^ ) 


2^ ±1 


X(C0/H2) 


2x10-^ 


10-5 


Tex (K) 


110 


110 


M„,oi (Mq) 


>0.28 [D/3kpc]2 


>0.38 [D/3kpc]2 



''A reliable estimate of this parameter is not possible due to low S/N 
''Integrated flux of the parabohc fit to the ^^CO hne profile (see Fig.lO & § 3.4) 
'^Affected by strong insterstellar absorption 
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Table 3: Envelope parameters from SED modeling 



Input source parameters 




Teff(K) 


7000 




6400 


D{pc) 


3000 


Input grain parameters 




Composition Astronomical silicates 


amin(ywni) 


0.01 


amax(/uni) 


0.3 


n(a) 


oc 


gas-to-dust mass ratio 


5=200 


Derived parameters 




Inner shell 




-Rin,i(cm) 


6.3x10^^ 


-Rout,! (cm) 


1.6x10^6 


UK) 


500-105 


p{r) 


_2 

OC r 


M{Mq yr-i) 


3.2x10-5 


M(Mo) 


0.007 


Outer shell 




i?in,2(cm) 


1.6x10^6 


-Rout,2(cm) 


1.25x10^^ 


UK) 


105-50 


p{r) 


oc r 


M{Mq yr-i) 


2.0x10-3 


M(Mo) 


0.95 


Large-grain-component 




a 


> 400^m 


UK) 


150-20 


M{Mq) 


0.03-0.3 
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Fig. 1. — Adaptive optics images of IRAS 18276—1431 in the Kp, Lp, and Ms near-infrared 
bands. A logarithmic intensity scale, from 3o" to 95% of the intensity maximum in each 
image has been used (see Table 1). Contours are for levels at 5%, and from 10% to 100% of 
the maximum in each image in 10% steps. 
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Fig. 2. — i>e/i(^ Kp-band image of IRAS 18276—1431 showing the faint, elongated halo beyond 
the bipolar lobes. Intensity scale (shown in reverse gray scale) is the same as in Fig.l. 
Contour levels are 3(7x2^ for i=0 to 11 ((7=8xlO~^ mJy pix~^, Table 1). The two bright 
objects at offsets — 1'.'4,+1'.'8 (referred to as FSl in the text) and +2^4,0 are not physically 
associated with IRAS 18276— 1431, as indicated by polarization maps of Gledhill (2005). 
Right) A sharpened version of the image on the left panel to show the "searchlight beams" 
and "arcs" of IRAS 18276— 1431. The maximum (black) and minimun (gray) values of the 
displayed logarithmic intensity scale are set at 0.8 and 1.3, respectively. The cross indicates 
the point where the searchlight beams intersect each other. The center of the image (offset 
0,0) is the center of the halo (see text in §3.1). 



Fig. 3. — A false three-color composite of the Kp (blue), Lp (green), and Mg (red) AO images 
in a squared-root intensity scale. Dashed lines indicate the orientation of the searchlight 
beams. Note the well delineated lobe caps and the equatorial waist of dust obscuring the 
central star. 
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Fig. 4.— F606W and F814W WFPC2/HST images of IRAS 18276-1431 in a logarithmic 
intensity scale. Contour levels are 60" x 2* for i=0 to 8 (see Table 1 for values of a). The 
North and East directions are indicated by the C/.'S-sized arrows. 
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Fig. 5. — Radial profile through the lobes along PA=17° of the Kp (black), Lp (blue), and 
Ms (pink) images showing the dark waist and the lobe caps. 
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Fig. 6. — Radial profile of the i^^p-sharpened image (Fig. 2 Right) through the searchlight 
beams integrating azimuthally in the range PA=6-10° and PA=34-38° showing the arcs 
(vertical tick marks). Question marks indicate tentative arcs. 
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Fig. 7. — Spectral energy distribution (SED) of IRAS 18276—1431 (squares) and the bright 
field star FSl (circles; see Fig. 2). Data from this paper are plotted using filled symbols. Lit- 
erature data-points are represented by open squares (from the 2MASS and IRAS databases, 
Le Bertre et al. 1987, van der Veen et al. 1995, Meixner et al. 1999). Crossed out symbols 
represent photometric data obtained with large apertures most likely including flux from 
FSl and that overestimate the flux from IRAS 18276—1431. A sohd (dashed) fine has been 
used to delineate the SED of IRAS 18276-1431 (FSl) in the 0.6-3.8//m range. 
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Fig. 8. — Top) Color maps obtained with the Kp, Lp, and Mg AO images. The color scale 
(in a logarithmic strecht) is shown in the wedge at the top of the boxes. Circles indicate 
the position and size of the small, relatively blue patches used for studying the central 
illuminating source (§ 5). Bottom) F814W— F606W color map (grey-scale); the F814W image 
(contours) is superimposed. 



Fig. 9.— ^^CO 1-0 velocity-channel maps of IRAS 18276-1431 obtained with OVRO. The 
LSR velocities of the different channels are indicated in km s~^ units within the boxes. Lev- 
els are —Aa, —2a, 2a, Aa... ((T=13Jybeam~^). Positive (negative) contours arc plotted 
using solid (dashed) lines. Brightness to temperature units conversion factor is 1.19 K per 
Jybeam^^. The 2.6mm continuum map and the clean beam, FWHM~ 10'.'9x7'.'l, are rep- 
resented in the last panel. Levels of the continuum are —2a, 2a, and Scr ((T=4mJy beam~^). 
The coordinates of the center of the maps (marked with a cross) are RA=18^30™30!75, 
Dec=-14°28'56'.'28 (J2000.0) 
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Fig. 10. — Left) ^^CO 1-0 spectrum integrated over the nebula (solid line) and parabolic fit 
to the ^^CO profile (dashed line). ^^CO (Middle) and ^^CO (Right) intensity maps integrated 
over the line profile, i.e., within the LSR velocity range ~[40:80km s~^] and [53:59]km s~^, 
respectively. For ^^CO (^'^CO) levels are —2a, 2a,... spaced every 2(T=3.4Jykm s~^(—2a, 
— la, la ... spaced every lcr=0.6Jykm s~^.) 
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Fig. 11. — Lp—Kp (left) and M^—Lp {middle) color maps (contours) superimposed on the 
sharpened i^'p-band image (grey-scale), and the Lp-band image {right). In the right panel 
only the lowest (highest) contours of the Lp—Kp {M^—Lp) map are shown. The cross plotted 
is located at the point where the beams cross at the nebula center and is oriented along at 
PA=17°. Note the similar orientation of the outer and inner equatorial structure probed by 
the Lp—Kp and M^-Lp color maps, and the Lp-band images. 
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Fig. 12. — Surface brightness radial profiles through the center of the halo along the nebular 
major axis (PA=15°-25°; top) and the equator (PA=105°-115°; bottom) of the Kp (black), 
Lp (blue), and Ms (pink) images. Cuts are shifted by an arbitrary amount in the relative 
intensity axis for clarity. Thick (thin) lines are used for the NE (SW) lobe. Radial power 
laws cx and cx r'^ are shown for comparison (dotted- lines). 
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Fig. 13. — Top) Observed SED of IRAS 18276—1431 corrected by interstellar extinction 
(squares), SED predicted by CLOUDY for our best-fit model (thin-solid, black line), and 
blackbody spectra for Tdust=150K and Tdust=20K (dashed hues). A combined model spec- 
trum (best-fit CLOUDY model plus 150 K black-body) is shown with the thick-solid, blue 
line. Bottom) Spatial variation of the grain equilibrium temperature for the different grain- 
size bins used in the model (solid, black lines). The gas density (blue line), obtained by 
adopting a constant dust-to-gas mass ratio thoughtout the envelope (5=200), and the gas 
kinetic temperature (pink, dashed line) predicted by CLOUDY for the best-fit envelope 
parameters derived from SED modeling are also shown. 



Fig. 14. — Model of spiral density shells induced in a circumstellar outflow by orbital motion 
in a binary system. The gray-scale shows the projected dust density for an observer lying in 
the system's orbital plane - projected as a horizontal line running through the center. (The 
projection from most other angles has a spiral form.) To compare with the observations of 
IRAS 18276—1431, the dust distribution should be convolved with the illumination pattern, 
which is strongest along a cone centered on the system's angular momentum vertical 
line through the center. The details of this model are described separately (Morris et al., in 
preparation). 



